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Preface 


The  objective  of  this  research  was  to  develop  a  primarily  physics-based  methodology  capable  of 
predicting  the  conductive  thermal  ignition  of  energetic  materials  by  residual  fragments.  The 
intended  application  of  this  methodology  is  to  provide  input  for  survivability  and  lethality 
analyses.  Our  research  used  a  combination  of  ballistic  range  experimentation  and  computer 
simulations  to  develop  the  methodology  described  in  this  report. 
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1.  Introduction 


With  the  current  emphasis  on  faster,  lightly-armored  vehicles,  small  and  medium  caliber-sized 
threats  are  a  greater  concern  than  they  would  be  with  traditional,  heavily-armored  vehicles. 

Aside  from  the  mechanical  damage  expected  from  perforating  threats,  there  is  an  added  concern 
about  the  vulnerability  of  energetic  materials  such  as  propellants,  explosives,  and  fuel  (1).  The 
ignition  of  energetic  materials  by  a  single  projectile  has  the  potential  to  cause  a  greater  amount  of 
large-scale  damage  than  a  single  projectile  is  capable  of  alone. 

The  threat  of  thermal  ignition  by  fragments  occurs  from  the  heating  induced  when  a  fragment 
perforates  armor  and  is  plastically  deformed.  As  a  fragment  penetrates  an  armor  plate,  it  is 
locally  heated  to  high  temperatures  in  the  regions  where  the  majority  of  plastic  work  has 
occurred.  When  these  heated  residual  fragments  come  into  contact  with  the  energetic  material 
behind  the  armor,  they  can  conduct  enough  heat  to  trigger  the  material  to  ignite  (2)  and, 
subsequently,  cause  significant  damage  to  personnel  and  vehicles.  Energetic  materials  can  also 
be  ignited  by  fragments  through  shock  compression,  if  the  residual  fragments  have  high 
velocities.  However,  this  research  pertains  to  relatively  low  residual  velocities  (Vr)  and  focuses 
on  the  conductive  thermal  ignition  regime. 

Energetic  materials  react  in  a  number  of  ways  to  ignition,  which  are  (in  ascending  order  of 
reaction  speed  and  violence):  bum,  fast  bum,  deflagration,  explosion,  partial  detonation,  and 
detonation  (5).  The  energetic  material  addressed  in  this  research  was  M30A1  tank  gun 
propellant.  In  our  experiments,  we  only  observed  reactions  in  the  bum  to  fast  bum  regimes. 


2.  Approach 


Our  research  focused  on  a  simplified  threat-target  pair  that  is  representative  of  medium  caliber 
penetrators  or  high  explosive  (HE)  munitions  fragments  and  light-armor.  The  penetrators  we 
studied  were  4340  steel,  right  circular  cylinder  (RCC)  fragment  simulating  projectiles  (FSPs) 
with  masses  of  208  and  830  grains  (gr),  or  0.5  caliber  (cal)  and  20  mm,  respectively.  The  ESP 
masses  are  representative  of  fragments  generated  by  large  HE  fragmenting  munitions.  The  target 
material  was  titanium,  precisely;  Ti-6A1-4V  extra-low  interstitial  (ELI)  annealed  plates  with 
thicknesses  of  0.25  and  0.5  in. 
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2.1  Ballistic  Experiment  Approach 


Two  types  of  ballistic  experiments  were  carried  out  to  support  this  research:  ignition  threshold 
velocity  shots  and  Vr  shots.  We  conducted  each  of  the  two  types  of  experiments  for  two 
combinations  of  FSP  mass  and  titanium  plate  thickness,  830  gr  FSP  vs.  0.25  in.  thick  titanium 
and  208  gr  FSP  vs.  0.25  in.  thick  titanium. 

2.1.1  Ignition  Threshold  Velocity  Experiments 

The  first  set  of  experiments  was  designed  to  measure  the  ignition  threshold  velocity  for  a  given 
FSP.  The  ignition  threshold  velocity  is  the  FSP  striking  velocity  (Vs).  We  observed  nearly  all 
the  higher  velocities  will  cause  ignition,  while  the  lower  velocities  will  not.  We  shot  the  FSPs  at 
0°  obliquity  through  the  titanium  plates.  The  titanium  plates  covered  the  only  open  side  of 
wooden  boxes  containing  M30A1  gun  propellant.  The  boxes  were  lined  with  fire-resistant 
gypsum  board  and  had  internal  dimensions  of  18x12x12  in.  A  0.5  in.  thick  square  of  Nomex  felt 
was  placed  directly  behind  the  titanium  plate  to  slow  the  residual  FSP  without  further  deforming 
it.  Slowing  the  FSP  after  perforation  helped  ensure  minimal,  if  any,  heating  was  done  in  the 
propellant  by  shock  deformation.  The  M30A1  propellant  grains  were  loosely  packed  in  the  box 
behind  the  titanium  plate  and  Nomex  felt.  Pieces  of  Celotex  and  high  density  foam,  usually 
totaling  6  in.  thick,  were  used  to  fill  the  remaining  volume  in  the  back  of  the  boxes.  The  FSP  Vs 
were  measured  with  radar,  and  a  video  camera  recorded  each  impact  and  subsequent  propellant 
reaction,  if  one  occurred.  Figure  1  shows  a  schematic  of  the  ignition  shot  experimental  setup  and 
figures  2a  and  2b  show  images  of  the  setup. 


Wooden  Box  Frame 


Loosely  Packed 
1  Propellant 
Brains 


Titanium 

Target 


Plate  Nomex 


Cellotex  and  High 
Density  Foam 


Ballistic 

Felt 


Figure  1 .  A  schematic  of  the  ignition  shot  experimental  setup. 
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Figure  2a.  An  image  of  a  box  with  its  top  removed. 


Figure  2b.  An  image  of  a  box  ready  to  be  shot. 


NOTE:  The  FSP  is  shot  through  the  hole  in  the  grey  steel  brace  on  the  front  of  the 

box.  The  brace  is  used  to  keep  the  box  in  place  in  the  event  of  a  propellant  reaction. 
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2.1.2  Residual  Velocity  Experiments 

Once  the  ignition  threshold  impact  velocities  were  determined,  the  resulting  Vr  of  the  fragments 
needed  to  be  established.  For  the  Vr  experiments,  we  shot  the  FSPs  through  the  titanium  plates 
and  measured  the  Vr  using  flash  X-ray  radiographs.  Each  fragment-plate  combination  was  shot 
at  its  observed  ignition  threshold  velocity.  The  plates  were  free-standing  without  the  boxes  from 
the  ignition  shots.  These  Vr  data  were  used  to  assist  in  validating  an  Eulerian  shock-physics 
code  (CTH)  computer  model  discussed  in  the  next  section. 

2.2  Computer  Simulation  Approach 

Due  to  time  restraints,  the  computer  simulations  focused  on  the  830  gr  FSP  vs.  0.25  in.  thick 
titanium  combination.  This  combination  was  chosen  because  it  provided  the  most  consistent 
results  and  presented  an  apparently  distinct  ignition  threshold  Vs. 

In  order  to  accurately  model  the  entire  armor  perforation  and  thermal  ignition  event,  we  used  two 
computer  codes,  CTH  and  a  Lagrangian  hydrocode  (LS-DYNA).  We  chose  CTH  to  model  the 
armor  perforation  event.  The  CTH  requires  the  use  of  the  Courant  stability  criterion  (CSC), 
which  requires  small  time  steps  in  CTH  (~10'^  s  for  the  armor  perforation  simulation).  The  small 
time  steps  make  heat  conduction  calculations  using  CTH  unreasonable  over  the  +2  s  required  to 
perform  the  entire  simulation.  The  LS-DYNA  allows  for  suppression  of  the  CSC,  and 
subsequently,  allows  for  longer  time  steps.  We  chose  to  use  LS-DYNA  to  run  the  heat 
conduction  calculations,  because  the  CSC  is  irrelevant  to  heat  conduction.  The  output,  from 
LS-DYNA,  allowed  us  to  predict  ignition  of  the  propellant. 

The  process  of  coupling  the  data  between  CTH  and  LS-DYNA  was  a  non-trivial  step.  The 
fundamental  differences  in  the  two  codes  made  it  difficult  to  translate  the  data  between  them. 
However,  a  simple  coupling  solution  was  developed  by  Raftenberg  and  Meyer,  and  is  described 
in  detail  in  references  4  and  5. 

Measurements  from  both  experimental  sets  were  used  to  feed  the  CTH  model  inputs.  The  three- 
dimensional  (3D)  CTH  penetration  simulation  modeled  the  impact  of  the  steel  FSP  against  the 
titanium  plate  at  the  ignition  threshold  Vs.  The  main  CTH  outputs  were  the  residual  geometry  of 
the  steel  FSP  and  its  3D  temperature  profile.  The  CTH  calculated  the  temperature  increment  of 
each  cell  in  the  mesh  by  dividing  the  work  increment  by  the  specific  heat  in  each  time  step.  The 
heat  conduction  was  neglected  in  these  CTH  calculations  because  it  was  irrelevant  over  the 
200  ps  long  event.  Figure  3  shows  the  steel  and  titanium  regions  of  the  CTH  simulation  at 
200  ps  after  impact. 
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Figure  3.  CTH  output  for  the  830  gr  steel  FSP  vs.  0.25  in.  titanium  plate  at  200  ps 

after  impact.  The  steel  material  is  colored  grey,  while  the  titanium  material 
is  blue  (4). 

In  a  late  review  of  the  data,  we  diseovered  the  density  and  dilation  eonstants  for  the  equation  of 
state  input  into  the  CTH  model  were  for  aluminum  instead  of  titanium.  All  of  the  other  model 
inputs,  including  those  for  strength  and  fracture,  were  correct.  However,  because  aluminum  is 
less  dense  than  titanium,  it  is  safe  to  assume  that  the  target  material  represented  in  our  CTH 
model  provided  slightly  less  ballistic  resistance  to  the  FSP  projectile  than  would  titanium  with 
accurate  equation  of  state  inputs.  A  detailed  description  of  the  inputs  used  is  provided  in 
reference  4  and  5. 

The  CTH  output  data  were  read  in  an  intermediate  processing  step  using  EnSight  and 
HyperMesh,  a  visualization  and  meshing  software,  respectively.  This  combination  of  software 
translated  the  3D  geometry  and  temperature  profde  data  into  a  format  acceptable  by  LS-DYNA. 
Again,  reference  4  and  5  discuss  this  intermediate  step  between  CTH  and  LS-DYNA  at  length.. 
The  LS-DYNA  simulation  modeled  a  single  M30A1  propellant  grain  in  intimate  contact  with  the 
leading  edge  of  the  FSP.  The  leading  edge  is  the  hottest  region  of  the  residual  FSP,  because  it  is 
the  region  that  experiences  the  highest  plastic  strain.  However,  because  our  analysis  used  a 
simplified  geometry,  the  propellant  was  not  in  contact  with  the  nodes  of  the  residual  FSP  that 
CTH  determined  to  be  the  absolute  hottest.  The  initial  conditions  for  the  LS-DYNA  simulation 
are  reproduced  in  figures  4a-4c. 
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Figure  4a.  The  steel  material  is  colored  orange,  while  the 
M30A1  propellant  grain  is  colored  blue. 


Figure  4b.  Surface  temperature  contours  of  the  LS-DYNA  initial  conditions.  The 
propellant  grain  is  room  temperature  and  the  FSP  is  heated  around  the 
plastically  deformed  leading  edge. 


Figure  4c.  Initial  temperature  contours  in  the  mid-plane  of  the  LS-DYNA  simulation. 
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The  LS-DYNA  heat  conduction  simulation  calculated  the  3D  thermal  conduction  between  the 
FSP  and  propellant  grain.  Deformation  and  motion  were  neglected  in  the  LS-DYNA  model 
because  we  assumed  that  the  residual  FSP  comes  to  rest  before  any  heating  occurs  in  the 
propellant.  The  heat  from  the  hottest  region  of  the  FSP  flowed  into  the  propellant  grain  and  into 
the  cooler  regions  of  the  FSP. 


3.  Results 


3.1  Experimental  Ballistics  Results 

The  data  presented  below  in  Table  1  display  the  results  of  the  ignition  threshold  velocity  shots 
for  the  208  gr  FSP  vs.  0.5  in.  titanium  plate  combination. 


Table  1.  Results  from  the  ignition  threshold  velocity  experiments  for 
the  208  gr  FSP. 


Shot 

No. 

FSP  Mass 
(gr) 

Ti64  Target 
Thickness 
(in.) 

V. 

(m/s) 

Burn? 

1 

208 

0.50 

807 

No 

2 

208 

0.50 

810 

No 

3 

208 

0.50 

819 

No 

4 

208 

0.50 

836 

No 

5 

208 

0.50 

839 

No 

6 

208 

0.50 

847 

Yes 

7 

208 

0.50 

848 

Yes 

8 

208 

0.50 

856 

Yes 

9 

208 

0.50 

869 

No 

10 

208 

0.50 

871 

No 

11 

208 

0.50 

878 

Yes 

12 

208 

0.50 

892 

Yes 

13 

208 

0.50 

913 

Yes 

In  shot  numbers  1,  2,  3,  4,  7,  and  9,  we  observed  a  phenomenon  where  the  residual  steel  FSP  was 
encased  by  a  titanium  plug  or  pieces  of  titanium  spall.  The  titanium  covered  the  FSPs  with 
varying  degrees  of  completeness.  In  shot  numbers  1,  2,  and  3  the  entire  residual  FSPs,  excluding 
the  rear  faces,  were  covered  by  titanium.  In  these  cases,  the  hot  leading  edges  of  the  FSPs  were 
completely  obstructed  from  direct  contact  with  the  propellant  by  titanium.  Figure  5a  shows  the 
residual  FSP  and  titanium  plug  from  shot  number  2,  which  is  representative  of  shot  numbers  1 
and  3  as  well.  In  shot  numbers  4,  7,  and  9,  a  titanium  spall  ring  covered  only  a  portion  of  the 
FSP,  leaving  the  front  face,  rear  face,  and  some  of  the  hot  leading  edge  exposed.  Figure  5b 
shows  shot  number  7,  and  indicates  where  the  heated  leading  edge  would  be  exposed  to  direct 
contact  with  propellant. 


7 


Figure  5  a.  Residual  FSP  from  shot  number  2. 

NOTE:  The  leading  edge  of  the  FSP  is  completely  covered  by  the  titanium  plug.  This 
residual  formation  is  also  representative  of  shot  numbers  1  and  3. 


Figure  5b.  Residual  FSP  from  shot  number  7. 

NOTE:  The  Ti  plug  and  steel  FSP  are  charred  black  because  this  residual  was  involved  in  a 

propellant  bum  reaction.  The  red  box  highlights  the  exposed  section  of  the  FSP’s  plastically 
deformed  leading  edge.  The  Ti  spall  ring  completely  covers  the  rest  of  the  leading  edge. 
This  residual  formation  is  also  representative  of  shot  numbers  4  and  9. 

Based  on  these  qualitative  observations  and  the  ignition  threshold  veloeity  data  in  table  1,  we  did 
not  settle  on  a  single  approximate  ignition  threshold  veloeity.  However,  the  data  seem  to  suggest 
that  it  lies  between  840  and  880  m/s. 

Figure  6  shows  the  results  from  the  Vr  shots.  Beeause  we  did  not  set  a  distinet  ignition  threshold 
veloeity  for  this  FSP-target  thiekness  eombination,  we  did  not  extrapolate  a  eorresponding 
distinet  Vr. 
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Striking  Velocity  vs.  Residual  Velocity 

208  gr.  vs.  .5"  Titanium 
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striking  Velocity  (m/s) 


Figure  6.  A  graph  of  data  from  the  Vj-  experiments  displays  Vs  vs.  the  corresponding  Vj-. 

The  inconsistencies  and  strange  plugging  phenomenon  observed  in  the  208  gr  FSP  set  of 
experiments  contributed  to  our  decision  to  focus  on  modeling  the  830  gr  FSP  set  of  experiments 
in  the  computer  simulation. 

The  data  presented  below  in  table  2  display  the  results  of  the  ignition  threshold  velocity  shots  for 
the  830  gr  FSP  vs.  0.25  in.  titanium  plate  combination. 


Table  2.  Results  from  the  ignition  threshold  velocity  experiments  for  the  830  gr  FSP. 


Shot 

Number 

FSP  Mass 

(gr) 

Ti64  Target 
Thickness 
(in.) 

V 

(m/s) 

Burn? 

1 

830 

0.25 

413 

No 

2 

830 

0.25 

424 

No 

3 

830 

0.25 

429 

No 

4 

830 

0.25 

429 

No 

5 

830 

0.25 

431 

Yes 

6 

830 

0.25 

432 

Yes 

7 

830 

0.25 

433 

No 

8 

830 

0.25 

434 

No 

9 

830 

0.25 

464 

Yes 

10 

830 

0.25 

482 

Yes 

11 

830 

0.25 

512 

Yes 
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The  perforating  830  gr  FSPs  also  indueed  plugging  in  the  0.25  in.  thiek  titanium,  but  the  plugs 
were  not  large  enough  to  eover  the  FSP,  and  they  were  always  clearly  separated  from  the 
residual  steel  fragments. 

In  shot  number  3,  smoke  was  observed  emitting  from  the  box  for  +15  s  before  it  eventually 
stopped  and  was  deemed  a  “No  Bum.”  Shot  5  was  also  notable,  in  that  it  took  a  longer  time  than 
most  shots  to  initiate  a  bum,  ~8  s  after  impact,  compared  to  an  ignition  time  of  <3  s  for  most 
other  shots.  Based  on  these  observations  and  the  data  in  the  table  2,  we  settled  on  an  ignition 
threshold  velocity  of  433  m/s  for  the  830  gr  vs.  0.25  in.  titanium  plate. 

Figure  7  shows  the  results  from  the  Vr  shots.  These  data  were  used  to  extrapolate  an 
experimental  Vr  for  our  ignition  threshold  Vs  of  433  m/s.  The  solid  green  line  indicates  where 
ignition  threshold  velocity  of  433  m/s  intersects  the  linear  fit  to  give  an  expected  Vr  of  ~205  m/s. 
This  data  was  used  to  help  confirm  the  validity  of  the  CTH  model. 


Striking  Velocity  vs.  Residual  Velocity 

830  gr.  vs.  .25"  Titanium 


Striking  Velocity  (m/s) 


Figure  7.  A  graph  of  data  from  the  Vf  experiments  displays  Vs  vs.  the  corresponding  Vr.  For  the  ignition  threshold 
Vs  of  433  m/s,  a  Vr  of  205  m/s  is  extrapolated  from  the  linear  fit  line. 
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3.2  Computer  Simulations  Results 
3.2.1  CTH 

The  results  reported  in  this  section  apply  only  to  the  830  gr  FSP  vs.  0.25  in.  thick  titanium  plate 
combination. 

The  CTH  calculated  a  Vr  of  215  m/s.  This  was  very  close  to  the  experimental  Vr  results, 

205  m/s,  reported  in  figure  7.  Figure  8  shows  the  X  dimension  velocity  of  the  FSP  as  calculated 
by  the  CTH.  The  X  dimension  is  positive  down  range;  normal  from  the  target.  The  Vr  is  taken 
as  the  velocity  of  the  FSP  at  200  ps  when  it  completely  exits  the  titanium  plate. 


Figure  8.  V,  of  the  830  gr  FSP  in  CTH  is  measured  at  200  ps,  when  the  FSP  exits  the 
titanium  plate. 

NOTE;  The  oscillations  observed  in  the  velocity  graph  are  due  to  shock  waves  reverberating 
in  the  FSP. 

The  residual  geometry,  of  the  FSP  produced  by  the  CTH,  showed  deformation  in  the  FSP’s 
leading  edge  that  was  consistent  with  the  deformation  geometry  observed  in  the  live 
experiments.  A  sample  comparison  between  a  real  FSP  and  the  CTH  geometry  is  shown  in 
Figure  9a  and  9b. 
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Figure  9a.  Photo  of  a  real  initial  FSP  (left)  and  real  residual  FSP  (right). 

NOTE:  The  residual  in  this  photo  was  recovered  from  the  ignition  threshold 

velocity  Shot  Number  7.  It  had  a  striking  velocity  equal  to  the  ignition 
threshold  velocity  of  433  m/s.  The  scale  in  the  photo  is  in  cm. 


Figure  9b.  An  image  of  the  residual  FSP  geometry 
produced  from  the  CTH  results. 

NOTE:  The  scale  of  the  image  frame  is  4  cm  vertically 
by  3  cm  horizontally.  The  small  stabilizing  fin 
on  the  rear  of  the  real  FSPs  was  not  included  in 
the  CTH  model. 

The  average  measured  bulge  diameter  of  the  residual  FSPs  shot  at  or  near  the  ignition  threshold 
velocity  of  433  m/s  was  ~2.25  cm,  a  12.5%  increase.  The  bulge  diameter  of  the  residual  FSP 
generated  by  CTH  was  2.17  cm,  an  8.5%  increase.  The  initial  FSP  diameter,  for  both  the  real 
and  simulated,  was  2.00  cm. 
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Figure  10  shows  the  temperature  fields  ealeulated  by  CTH  for  a  residual  FSP  with  Vs  of  433  m/s. 


Fringe  Levels 
5.0000402 
4.00004  02. 
4.0000402 
4.^1000402 

4.2000402 
4.0000402. 
3.0000402 
3.6000402 
3.il00ei0Z 

3.2000402 
3.0000402 


Figure  10.  Temperatures  within  the  midplane  of  the  residual  FSP  as 
calculated  by  CTH  at  t  =  200  ps  {4). 

NOTE:  The  hottest  region  is  the  leading  edge  of  the  RCC,  which  underwent 

the  most  plastic  deformation,  while  the  rest  of  the  FSP  is  relatively  cool. 
These  results  are  applied  as  initial  conditions  in  the  LS-DYNA  analysis. 
(The  fringe  levels  are  kelvin.) 
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3.2.2  LS-DYNA 


The  CTH  results  shown  in  figure  10  are  the  initial  eonditions  for  the  residual  FSP  in  LS-DYNA 
which  are  also  shown  in  figure  4  along  with  the  adjacent  propellant  grain. 

Figure  1 1  shows  the  initial  temperature  field  on  the  face  of  the  propellant  grain  in  direct  contact 
with  the  steel  residual,  face  A  (x,  z  plane).  Face  B  is  the  plane  of  propellant  material  recessed  one 
element  dimension,  0.635  mm,  into  the  grain  and  along  the  axis,  from  face^.  Nodes  lA 
through  6A  are  the  six  hottest  nodes  in  face  A  and  are  shared  by  the  steel  residual.  Nodes  IB 
through  6B  are  their  counterparts  with  the  same  x  and  z  coordinates  but  on  face  B. 


Fringe  Levels 


5.000e+02 


4.800e+02 


4.600e+02 


A.AUUe+UZ 


4.200e+02 


AmUe+UZ 


3.800e+02 


3.688e+82 


3.4BBe+B2 


3.2flfle+fl2 


3.BBBe+B2 


Figure  1 1 .  The  initial  temperature  field  in  the  M30A1  grain  and  the  six-nodes  and  two-faces  (A  and  B)  for  which 
temperature  histories  are  recorded  in  figure  12(4). 


Heat  from  the  highly  deformed  region  of  the  FSP  flows  faster  into  the  unheated  regions  of  the 
FSP  than  into  the  adjacent  propellant  grain.  This  is  due  to  the  higher  thermal  conductance  in  the 
steel,  38.1 1  W/m-K  (6),  compared  to  a  0.418  W/m-K  (7)  conductance  in  the  M30A1.  This  effect 
can  be  observed  in  figure  12.  Ignition  occurs  if  enough  heat  is  drawn  into  the  propellant  grain  to 
reach  an  ignition  threshold  temperature  before  the  unheated  region  of  the  FSP  conducts  it  away. 
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Figure  12.  Temperature  contours  in  the  midplane  of  residual  FSP  and  propellant  grain  at 
1  and  2  s  after  initial  contact  (4).  The  fringe  levels  are  kelvin. 
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Figure  13  shows  a  representative  temperature  history  of  the  hottest  regions  in  the  M30A1 
propellant.  The  location  of  nodes  6A  and  6B  are  defined  in  figure  11.  Node  6A  begins  at  524  K. 
At  2.4  s,  temperatures  of  nodes  6A  and  6B  have  converged  to  328  K  and  308  K,  respectively, 
suggesting  that  no  further  heating  occurs  in  the  propellant  after  2.4  s. 
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Figure  13.  Temperature  (K)  vs.  time  (s)  at  two  representative  nodes  identified  in  figure  11,  one  on  face 
A  and  one  on  face  B  (4),  and  ignition  times  reported  in  Rocchio  and  Wires  (2). 


The  data  points  in  figure  13  labeled  “R&W  Lower”  and  “R&W  Upper”  refer  to  M30  ignition 
data  reported  by  Rocchio  and  Wires  (R&W)  (2).  Rocchio  and  Wires  performed  temperature 
controlled  experiments  on  single  grains  of  M30.  For  a  range  of  temperatures  between  483  K  and 
533  K,  the  time  required  for  ignition  was  measured.  Their  results,  after  converting  from  Celsius 
to  kelvin,  are  reproduced  in  table  3.  In  our  experiments,  most  ignitions  occurred  2-3  s  after 
perforation  of  the  titanium.  However,  some  ignitions  occurred  at  upwards  of  8  s. 
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Table  3.  Results  from  the  R&W  temperature  controlled  ignition 
experiments  for  single  grains  of  M30  propellant  grains. 


Temperature  (K) 

Slower  (s) 

tupper  (s) 

533 

0.260 

0.378 

523 

0.574 

0.621 

513 

1.538 

2.329 

503 

3.422 

4.111 

493 

6.257 

7.045 

483 

9.400 

18.256 

To  quote  from  R&W,  “the  high  limit  was  that  time  above  which  one  could  always  expect 
sustained  combustion  and  the  low  limit  was  that  time  below  which  sustained  combustion  did  not 
occur”  (p.  14).  These  data  are  the  basis  for  an  initiation  criterion  to  apply  to  the  LS-DYNA 
temperature  solution.  In  figure  13,  since  all  points  on  curves  6 A  through  6B  lie  below  the  R&W 
data  points,  one  can  conclude  that  the  M30A1  grain  that  contacts  the  steel  residual  in  the  location 
and  configuration  depicted  in  figures  4  and  12  will  not  undergo  sustained  combustion. 


4.  Discussion 


Temperatures  in  the  M30A1  were  expected  to  come  close  to  the  ignition  data  in  R&W  (2) 
because  experiments  conducted  here  indicated  that  a  Vs  of  433  m/s  resulted  in  threshold  ignition 
of  the  propellant.  The  results  reported  here  suggest  that  the  methodology  under-predicts  the 
temperatures  generated  in  the  M30A1  by  the  hot  steel  residual. 

Several  uncertainties  in  our  analysis  may  have  contributed  to  the  model’s  under-prediction  of  the 
propellant  temperatures.  The  CTH  model  used  density  and  dilation  constants  for  the  target  plate 
that  were  representative  of  aluminum  instead  of  titanium.  This  may  have  contributed  to  the 
under  prediction  of  plastic  deformation  in  the  FSP  calculated  by  CTH  (8.5%  increase  in  diameter 
vs.  12.5%  for  experimental,  see  figure  9).  It  is  possible  that  with  the  proper  titanium  constants, 
CTH  would  have  calculated  more  deformation  and  higher  temperatures  in  the  residual  FSP.  This 
would  translate  to  higher  temperatures  calculated  by  LS-DYNA  within  the  propellant.  Our 
analysis  also  used  a  simplified  geometry  in  LS-DYNA  where  the  propellant  grain  was  in  contact 
with  the  heated  region  of  the  residual  FSP,  but  not  with  the  absolute  hottest  node  in  the  FSP’s 
heated  leading  edge.  Additionally,  the  temperature  computing  algorithms  in  CTH  have  a  limited 
predictive  capability;  therefore,  we  were  unable  to  experimentally  verify  the  temperature  in  the 
residual  FSP. 
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Despite  the  uncertainties  in  our  analysis,  the  methodology  is  promising  and  has  sufficient 
flexibility  to  model  other  materials  and  provide  ignition  predictions,  given  accurate  inputs  and 
ignition  temperature  data. 


5.  Conclusions 


The  results  of  two  sets  of  experiments  are  presented  in  this  report.  These  experiments  measured 
the  ignition  threshold  Vs  and  Vr  for  two  different  steel  RCC  FSP  masses,  830  gr  and  208  gr, 
targeting  two  thicknesses  of  titanium  plate  armor,  0.25  in.  and  0.5  in.,  respectively.  We  used 
these  experiments  to  develop  a  computer  simulation  method  for  modeling  and  predicting  the 
conductive  thermal  ignition  of  an  energetic  material  by  the  residual  steel  penetrators.  Our 
method  employed  two  computer  codes  in  series,  CTH  and  LS-DYNA.  The  codes  were  coupled 
using  a  intermediate  process  to  input  the  temperature  field  results  produced  by  CTH  into 
LS-DYNA  for  the  heat  conduction  analysis.  Using  this  method,  we  were  able  to  model  the  entire 
event  from  the  armor  perforation  to  the  propellant  ignition. 

The  method  was  used  to  model  the  perforation  of  a  0.25  in.  thick  titanium  plate  by  an  830  gr  FSP 
with  a  Vs  of  433  m/s  and  the  subsequent  contact  of  an  M30A1  propellant  grain  with  the  hot  steel 
residual.  Heat  conduction  was  neglected  in  the  steel-titanium  perforation  analysis,  performed 
with  CTH.  Deformation  and  motion  were  neglected  in  the  post-processing  heat  conduction 
analysis  performed  with  LS-DYNA.  The  temperature  field,  as  a  function  of  time  within  the 
M30A1  grain,  was  compared  with  data  in  R&W  (2)  for  ignition  time  as  a  function  of 
temperature.  The  temperature  field  in  the  M30A1  grain  never  reached  the  ignition  temperature 
data  in  R&W  (2),  indicating  a  prediction  of  no  ignition. 
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Abbreviations  and  Acronyms 


3D 

three-dimensional 

ARL 

U.S.  Army  Research  Laboratory 

cal 

caliber 

CSC 

Courant  Stability  Criterion 

CTH 

an  Eulerian  shock-physics  code 

ELI 

extra  low  interstitial 

ESP 

fragment  simulating  projectile 

gr 

grain 

HE 

high  explosive 

LS-DYNA 

a  Lagrangian  hydrocode 

R&W 

Rocchio  and  Wires 

RCC 

right  circular  cylinder 

SLAD 

Survivability/Lethality  Analysis  Directorate 

Vr 

residual  velocity 

Vs 

striking  velocity 
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A  KEINTZ 

WSMRNM  88002-5502 
USARL 

AMSRD  ARE  SL  E 
R  FLORES 

WSMRNM  88002-5513 

US  ARMY  EVALUATION  CTR 
CSTE  AEC  SVE 
R  LAUGHMAN 
4120  SUSQUEHANNA  AVE 
APG  MD  21005-3013 
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No.  of 
Copies 

42  HCs 
4  CDs 
1  elec 


Organization 

US  ARMY  RESEARCH  LABORATORY 
AMSRD  ARE  WM  TD 
M  RAFTENBERG  (5  HCs,  1  CD) 

H  MEYER  (5  HCs,  1  CD) 

T  BJERKE 
AMSRD  ARE  SL  BE 
L  ROACH 
E  FIORAVANTE 
M  MAHAFFEY 
J  ABELL 

K  BATES  (5  HCs,  2  CDs,  1  elec) 
AMSRD  ARE  SL  BM 
D  FARENWALD 
M  PERRY 
R  SAUCIER 
A  DIETRICH 
A  PRAKASH 
A  MIKHAIL 
T  MALLORY 
G  BRADLEY 
C  BARKER 
D  LYNCH 
G  MANNIX 
AMSRD  ARE  SL  BS 
W  WINNER 
D  HISLEY 
E.O.  DAVISSON 
M  GILLICH 
J  AUTEN 
M  KUNKEL 
AMSRD  ARE  SL  BD 
R  GROTE 
F  MARSH 
P  DOMINICK 
LMOSS 

AMSRD  ARE  HR  MB 
P  FEDELE 
APG  MD  21005-5066 


No.  of 

Copies  Organization 

14  HC  DIR  US  ARE 

AMSRD  ARE  SL 
J  BEILFUSS 
J  FRANZ 
P  TANENBAUM 
AMSRD  ARE  SL  BM 
M  PERRY 

AMSRD  ARE  SL  BB 
S  SNEAD 

AMSRD  ARE  SL  BD 
R  GROTE 

AMSRD  ARE  SL  BE 
L  ROACH 

AMSRD  ARE  SL  BM 
D  FARENWALD 
AMSRD  ARE  SL  BS 
W  WINNER  (4  HCS) 

AMSRD  ARE  SL  EC 
J  FEENEY 

I  HC  ASST  SECY  ARMY 

ACQSTN  LOGISTICS  &  TECH 
SAAL  ZP  RM  2E66I 
103  ARMY  PENTAGON 
WASHINGTON  DC  20310-0103 

I  HC  ASST  SECY  ARMY 

ACQSTN  LOGISTICS  &  TECH 
SAAL  ZS  RM  3E448 
103  ARMY  PENTAGON 
WASHINGTON  DC  20310-0103 

TOTAL:  82  (2  Electronic,  12  CDs,  and  68  HCs) 
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